Interplay of fixed points in scalar models 
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We performed the renormalization group analysis of scalar models exhibiting spontaneous sym- 
metry breaking. It is shown that an infrared fixed point appears inducing a dynamical scale, that 
can be identified with the correlation length of the model. This enables one to identify the type of 
the phase transition which shows similarity to the one appearing in the crossover scale. 
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/. Introduction. — The renormalization group (RG) 
method [1] is one of the best candidate to take into ac- 
count the degrees of freedom of a quantum field theo- 
retical model systematically. The method enables us to 
identify the fixed points of the models and the scaling 
of the couplings in their vicinity which can give us the 
critical exponents of the corresponding fixed point [2- 
4]. The RG method is usually tested in 3-dimensional 
(3D) 0{N) scalar models, where a trivial Gaussian and 
a non-trivial Wilson-Fisher (WF) fixed point exist. The 
calculation of the critical exponents of the latter fixed 
point plays the test ground of all inventions in the RG 
method [2, 5-8]. The scaling of the correlation length ^ 
defines the exponent v as 



(1) 



with the reduced temperature t. There is a second order, 
Ising type phase transition in the 3D 0(1) model with 
exponent v — 0.63. The Kosterlitz-Thouless (KT) or in- 
finite order phase transition [9, 10] is widely investigated 
too, furthermore it gives a great challenge to recover the 
scaling of ^ according to 



log ^ oc t 
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A typical example for KT transition is presented by the 
sine-Gordon (SG) model [4, 11-13] in two-dimensional 
(2D) Euclidean space, which belongs to the same univer- 
sality class as the 2D Coulomb gas and the 2D XY model. 
Furthermore the SG model and its generalizations with 
compact variable have been thoroughly investigated in 
the framework of integrable field theory [14]. 

The fixed points and the typical scaling of ^ in their 
vicinity as in Eq. (1) for the 3D 0{N) and in Eq. (2) for 
the 2D SG models are located in the crossover region, i.e. 
between the ultraviolet (UV) and infrared (IR) regimes. 
It has been argued [15] that there exists an IR fixed point 
in the 3D 0{N) model at a finite momentum scale, which 
can be uncovered by a genuine rescaling of the couplings 
around the singularity of the RG evolution, i.e. where the 
flows stop. Recently it has been obtained that there ex- 
ists a similar IR fixed point in the SG model, too [4]. We 
note that both IR fixed points belong to the spontaneous 
symmetry breaking phase of the models. 



We show in this letter that this IR fixed point is ac- 
companied by a dynamical length scale. It is defined by 
the degeneracy of the RG evolution, and can be identi- 
fied with the reciprocal of the correlation length. This 
provides us a new method to determine the exponent z/ 
beyond the crossover scaling regime in the vicinity of the 
IR fixed point. The IR fixed point can also account for 
the scaling of ^ and gives similar results that are ob- 
tained around the crossover fixed point of the model. 
The method is also applicable when there is no crossover 
fixed point, e.g. in the bi-layer SG model [16], showing a 
greater flexibility and a more fundamental nature of our 
method presented here. 

//. Evolution equation. — The Wetterich RG equation 
for the effective action is [17] 



T, = 
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where ■ = kdk, ' — d/d^p and the trace Tr denotes the 
integration over all momenta. Eq. (3) has been solved 
over the functional subspace defined by the ansatz 

T 
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with the 0(1) symmetric or the periodic potential Vk 
and the wave- function renormalization Z^, which con- 
stant part is z=Zk{^p=Q). The polynomial suppression 
has the form R = p^{k'^/p'^f, with 6 > 1. Eq. (3) leads 
to the evolution equations for the couplings [1, 5]. The 
loop integral appearing in the RG equations should be 
performed numerically. The scale k covers the momen- 
tum interval from the UV cutoff A to zero. Typically we 
set A = 1. If we introduce k — min(2:p^ + i?), the RG 
evolution becomes singular at k = kf when 



0, 
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where fc^ = bk^[z/{b - l)]l-l/^ The solution of this 
equation defines the scale at which the potential becomes 
degenerate. 

//. Ising type transition. — For the 3D 0(1) model the 
potential in Eq. (4) has the form 
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where p = </'^/2 and the last cquahty holds up to an irrel- 
evant additive constant. Higher precision can be achieved 
by the latter form of the potential, where the number 
(n — 1) couplings and the additional running parameter 
K figure. One can take into account the evolution of the 
wave function renormalization with similar ansatz for Zk 
as 



(7) 



where a similar p dependent ansatz can be formulated. 
We also use the normalized couplings defined 
x/lc', where x can be g2i, k, Z2i, i = 1, . . . n. Even bet- 
ter results can be obtained for the critical exponents 
without Taylor expanding the potential and the wave- 
function renormalization, but at this present state it is 
not required. The main point here is to demonstrate the 
equality of the exponent v taken from the vicinity of the 
WF fixed point and from the IR fixed point. In the 3D 
0{N) model the exponent v is identified as the negative 
reciprocal of the single negative eigenvalue of the matrix 
coming from the linearization of the evolution equations 
around the WF fixed point. In the lowest order of the 
gradient expansion, i.e. in the local potential approxima- 
tion (LPA) , when Zk = 1, the exponent i/ at the WF fixed 
point is 1/ ~ 0.53 {v w 0.64) for n=2 (n=4), respectively 
[6] . We always investigated the scheme-dependence of the 
results by choosing different values of b. The LPA fiows 
run into singularity at a finite scale kf due to the strong 
truncation of the flow equations. The reduced tempera- 
ture is identified as t ~ {ka — Kc) [7], i.e. the deviation of 
the UV initial coupling k to its critical value Kc, which is 
taken on the separatrix between the WF and the Gaus- 
sian fixed points. One can express the coupling k by the 
couplings g2n, e.g. for n=2 one has k ~ 92/94- The re- 
sults of the flow show that the parameter k blows up, i.e. 
At ^ cx) at the same value of the scale kf [15, 18], which 
characterizes the degenerate potential with its minimum 
in the infinity, and the building up of a global conden- 
sate [19]. We identify the inverse of A;/ as the correlation 
length ^. This enables us to read off the exponent u 
around the IR fixed point. The LPA results can be seen 
in the inset of Fig. 1 for n = 2 and n = 4 couplings, which 
shows the coincidence of the exponents calculated from 
the data around the WF and the IR points, and demon- 
strates the possibility for finding this quantity in the IR. 

By taking into account more and more couplings and 
including the evolution of Z^ one can push the singu- 
larity into k=0 giving a degenerate effective potential in 
the IR. Fig. 1 shows the phase diagram of the D3 0(1) 
model projected onto the normalized couplings 92 and §4, 
with the parameter choice 6 = 2. Similar phase structure 
can be obtained for any value of b. We numerically de- 
termined the exponent u for n = 8 couplings from data 




FIG. 1: The phase structure of the D3 0(1) model is pre- 
sented, where the evolution of Zk is taken into account. The 

upper cud of the thick curve denotes the WF fixed point. The 
trajectories tending to the right (left) correspond to the sym- 
metric (symmetry broken) phase, respectively. The scaling of 
the correlation length as the function of the reduced tempera- 
ture is plotted in the inset. The points are obtained from the 
scaling around the IR fixed point, while the solid lines rep- 
resent the scaling around the WF fixed point, i.e. v = 0.53 
(LPA), u = 0.64 (LPA), and u = 0.62 (including the flow of 
Zk) ioT n — 2, n = 4, and n = 8, respectively. The curves are 
shifted for better visibility. The circle and square correspond 
to 6 = 2, 5, respectively. 



around the WF and the IR points, which, similarly to the 
LPA results, shows high coincidence, as is demonstrated 
in the inset of Fig. 1. The numerical results also show 
that the wavefunction renormalization constant z blows 
up in the vicinity of the degeneracy as the function of k, 
as is demonstrated in Fig. 2. The blowup of the fiow of 




FIG. 2: The scaling of the z for n = 8 and b = 2 for several 
UV vahies of (74. The dotted line corresponds to the flow 
along the separatrix. In the deep IR regime the flows blow up 
in the broken symmetric phase, while they run into constant 
values in the symmetric one. 
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z occurs at smaller and smaller values of fcy as the UV 
initial value ka approaches its critical value. Thus the 
parameters z and k show up quite similar singular be- 
haviour. The blowup of z results in the stop of the flows 
of all the couplings at fc/ in the broken symmetric phase, 
while in the symmetric phase z goes to a constant value, 
giving LPA evolution in the IR. The flow of z denoted by 
dotted line in Fig. 2 correspond to the flow along the sep- 
aratrix in Fig. 1. Its slope gives the known, tiny anoma- 
lous dimension r] = — rflogz/dlogfc « 0.05 belonging to 
the WF fixed point. The flows in the vicinity of the WF 
fixed point pick up the effects coming from the fluctua- 
tions affected by the WF fixed point and bring them to 
the IR fixed point. There the anomalous dimension is ex- 
tremely high (and scheme dependent), but the exponent 
V must not change, since it characterizes the condensate 
of the broken phase. Thus, well beyond the WF fixed 
point scaling regime, in the vicinity of the IR fixed point 
in the deep IR one can recover all the information on the 
system. 

Periodic model. — The 2D SG model is defined via 
the potential 




Vk = Uk cos <p. 
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FIG. 3: Phase diagram of the SG model, with 6 = 5. The 
dashed (soUd) lines represent the trajectories belonging to 
the (broken) symmetric phase, respectively. The wide line 
denotes the separatrix between the phases. The inset shows 
the scaling of the correlation length as the function of the 
reduced temperature t. The curves are shifted for better vis- 
ibility. The lower (upper) set of lines corresponds to the IR 
(KT) fixed point. The triangle, circle and square correspond 
to 6 = 2, 5, 10, respectively. In the middle a straight line with 
the slope —1/2 is drawn to guide the eye. 



The higher harmonics of the SG model are neglected. 
They correspond to vortices with higher vorticity of the 
equivalent gas of topological excitations. It is known that 
only the fundamental mode plays a significant role in 
the determination of the thermodynamical properties of 
the model, while effects corresponding to higher vorticity 
are negligible. We also note that the fundamental mode 
can recover the phase structure of the model including 
the KT transition point [11] and the IR behaviour [4]. 
Furthermore the wavefunction renormalization constant 
z can account for the KT transition [4, 10], therefore we 
omit further terms in Z^- In the LPA approximation, 
the RG treatment of the SG model in the IR shows two 
phases separated by the Coleman point and the occur- 
ring the dynamical momentum scale where the cvohition 
of the normalized coupling u = u/k^ becomes marginal 
[12]. Identifying this scale as the reciprocal of the correla- 
tion length we obtain v ~ 1. The Coleman point becomes 
the KT transition point [4, 11] if evolves. Furthermore 
an additional IR fixed point turns up that can be trans- 
formed to the unique point at u = 1 and 1/z when 
the normalized coupling u is made of use [4]. Then any 
choice of h gives qualitatively similar phase diagrams. We 
plotted the case & = 5 in Fig. 3. The normalized coupling 
u tends to 1 for every value of b. It shows that the de- 
generate potential occurs in the IR limit of the broken 
symmetric phase independently of the RG scheme. This 
reflects the serious limitation of the LPA results. In the 
symmetric phase the evolution of z is negligible giving 
the same evolution as was obtained in LPA with the line 
of fixed points. 



One can easily show that the critical exponent 77,, char- 
acterizing the vortex-vortex correlation function [20] is 
?7„ = 1/4 independently of the parameter h [4]. However 
the anomalous dimension being characteristic for the di- 
vergence of the correlation function of the field variables 
gives ?7 = in the vicinity of the KT point. In the deep 
IR scaling region the situation changes significantly, there 
new scaling laws appear. Fig. 4 shows that around the 
KT point (at about k/\ ^ 10^^) z is practically con- 
stant, giving ?7 = 0, while in the IR region z diverges with 
the exponent ry = 26/(6 — 1), giving scheme-dependent 77 
values in the IR. It implies that one cannot avoid the 
evolution of z by including a phase factor z = k'^ m the 
LPA evolution equations. The scale k makes the evolu- 
tion of u and z qualitatively scheme-independent. The 
value of the wave-function renormalization z blows up 
at a certain minimal value of k in the broken symmetric 
phase, similarly to the 3D 0(1) model. Likewise, the IR 
fixed point of the broken symmetric phase is reached at 
A; = 0, but the IR scaling occurs already at some finite 
dynamical scale fc/. 

The initial UV value of za can be identified with the 
square of the temperature, thus its distance from the 
separatrix (for fixed ua) gives the reduced temperature 
t. Ast ^ Q the correlation length increases as in the case 
of the 3D 0(1) model. Our numerical results are shown 
in the inset of Fig. 3. There are two types of correlation 
lengths, one is defined as usual, namely at around the KT 
turning point of the coupling u. Another one is identified 
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FIG. 4: The evolution of z for different values of 6. The inset 
demonstrates that r\ = 2b/{b — 1) in the deep IR region. 



as ^ = in the neighbourhood of the IR fixed point. 
It can be seen from the inset of Fig. 3 that the scaling of 

^ shows an infinite order phase transition, for all schemes 
with the exponent « 0.5, as opposed to the LPA result 
i/= 1. 

Such a definition ^ in the IR is also powerful when we 
have no crossover scaling. In the case of the bi-layer SG 
(LSG) model [16] one can easily show that the model 
has no fixed point because of the evolution of the intcr- 
layer coupling. However, the detailed RG investigations 
performed by us have shown that there is an attractive 
IR fixed point of the LSG model with the appearing dy- 
namical scale kf as in the SG model, and gives the same 
KT type phase transition with exponent v « 0.57 [21], 
proving the infinite nature of the phase transition there. 

The scaling laws at the IR fixed point seem to coincide 
with those at the crossover. The IR fixed point is the 
signal of the spontaneous symmetry breaking in both the 
polynomial and the periodic models, while the KT and 
'WF points are crossover fixed points, although they are 
closely related to the IR one. The dependence of the cor- 
relation length ^ on the reduced temperature has already 
been picked up by the trajectory in the crossover regime 
and is carried by it to the IR fixed point. In that sense 
there is an interplay between the IR and crossover points, 
but the proper IR physics (including both the degener- 
acy and the scaling of the correlation length) can only 
be deduced from the information on the IR fixed points. 
This finding clearly demonstrates the fundamental prin- 
ciple of the RG method, namely the IR physics should 
recover all the information on the system. 
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